Grand Unification II: Hot Accretion and AGN Jets 

Renyue Cen^ 

Received ; accepted 



^Princeton University Observatory, Princeton University, Princeton, NJ 08544 
cen@ astro .princeton. edu 



- 2 - 



ABSTRACT 

We show that a quasi-spherical (QS) hot accretion flow is expected to operate 
in all SMBHs, with its rate being capped at rngs^max ~ 0.001 (M/10^ Mq) in units 
of the Eddington rate. It is then proposed that AGN jet power is proportional 
to the product of the hot accretion rate and a SMBH spin-dependent efficiency 
for energy extraction in the form of jets. Predictions from this model include 
(1) while radio jets should emerge from all SMBHs, the maximum jet power 
goes with SMBH mass approximately as 10^^-^(M/10® M0)^erg/s, (2) even 
although there are two separate underlying populations of radio-quiet (RQ) and 
radio-loud (RL) AGNs, any bimodality in the observed AGN radio loudness 
distribution is likely due to a selection effect, such as some imposed optical 
magnitude limits, (3) the RL fraction of quasars is expected to decrease with 
redshift in the cold dark matter model, (4) host galaxies of RQ and RL quasars 
should be drawn from the same underlying elliptical galaxy population, although 
RL quasars may have SMBHs that are somewhat more massive than their 
RQ counterparts and RL quasars may reside predominantly in core elliptical 
galaxies, (5) RL low-luminosity AGNs and LINERs may represent the long, 
declining "trailing" phase following the initial, more luminous AGN phase, 
(6) a broad anti-correlation between radio-loudness and disk accretion rate is 
expected, (7) RL AGNs may be expected to be more abundant in Type la 
supernovae than RQ AGNs, (8) among RL AGNs a correlation between radio 
power and clustering strength is predicted, and (9) RQ and RL AGNs should 
have, on average, similar IR-optical-UV properties. 

Subject headings: black hole physics - galaxies: active - galaxies: bulges - 
galaxies: interactions - radio continuum: galaxies - quasars: general 



-3- 



1. Introduction 

The past decades have witnessed a tremendous amount of effort devoted to 
understanding the observed complexities of AGNs. This pursuit of intrinsic simphcity 
has been remarkably fruitful. The most successful unification scheme hinges on a simple 
assumption that the observed optical/UV properties of an AGN depends only on the 
relative orientation of the observer's line of sight to the non-spherical configuration of the 
central engine, which is believed to be composed of a SMBH (Rees 1984), an optically thick 
and geometrically thin accretion disk, a surrounding hot corona, an obscuring torus of gas 
and dust (Antonucci 1993) and a diverse set of gas clouds moving at a variety of speeds 
at different locations (e.g., Elvis 2000). Unification proposals among RL AGNs, based on 
viewing orientation dependent Doppler effect of the relativistic (twin) radio jets, also enjoy 
great successes (Barthel 1989; Urry & Padovani 1995). 

The observational fact that most of, and perhaps all, non-active, bulged galaxies 
harbor SMBHs at their centers (e.g., Richstone et al. 1998) and the existence of a tight 
correlation between the observed central SMBH mass and bulge stellar mass (or velocity 
dispersion) (e.g., Magorrian et al. 1998; Fcrrarese & Merritt 2000; Gebhardt et al. 2000; 
Merritt & Ferrarese 2001; Tremaine et al. 2002) strongly suggest that growth of SMBHs 
and formation of galaxies, in particular, the galactic spheroid component, are intricately 
linked (e.g., Haehnelt, Natarajan, & Rees 1998). Motivated by these strong observational 
indications, we intend to embark on, through a series of focused papers, a coherent physical 
framework for coupled formation and evolution of galaxies and AGNs, in hopes of ultimately 
helping build a grand unification theory. Earlier attempts along this direction are typified 
by the pioneering work of Kauffmann & Haehnelt (2000), based on plausible but heuristic 
assumptions on coupled formation of galaxies and AGNs. In the first paper of this series 
(Gen 2006), we put forth a more physically provoking explanation for the observed SMBH 
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mass-bulge velocity dispersion relation, in the context of hierarchical structure formation 
model. This paper is the second of the series. We show a co-existence of two accretion (hot 
and cold) modes on SMBH scales and propose that the hot accretion, in conjunction with 
the spin of the SMBH, determines the AGN jet power. In §2 we present a theory on the 
dependence of hot accretion rate and the jet power on the mass of the SMBH. In §3 some 
observational tests and predictions are made, followed by conclusions in §4. 



We are interested in knowing the state of hot gas on its way to the SMBH. For ensuing 
derivations, we utilize the observed SMBH mass (M)-bulge (1-d) velocity dispersion (a) 
relation (Tremaine et al. 2002): 



Cen (2006), in a specific model, showed that this relation may not evolve with redshift. 
Assuming that the gas temperature T relates to a by kT/fj, = rUpa'^, where fj, = 0.58 is 
the molecular weight, rUp the proton mass and k the Boltzmann constant, we obtain (with 
specific heat 7 = 5/3) 



With Equations (1,2) and the adopted temperature- velocity dispersion relation we can 
further relate the gas density at r^, ub, to the (adiabatic) Bondi (1952) accretion rate, rriB, 
in Eddington units (Me), as a function of M: 



2. Hot Accretion and AGN Jet Power: Theory 



2.1. Basic Assumptions and McLximum Hot Accretion Rate 




(1) 




(2) 




(3) 
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where Me = 2.2(M/10^Mq) Mg/yr is the Eddington rate. Equation (3) provides a way to 
relate hot gas coohng time to the accretion rate for a fixed SMBH mass. In what follows all 
accretion rates are expressed in the Eddington rate. 

We distinguish a hot, quasi- spherical Bondi accretion mode of rate TfiQs from a cold 
disk accretion mode. While the cold disk accretion mode of rate, rhd = rhhc + rhc, consist of 
two components, one of rate rhhc from cooling of a hot component en-route from galactic 
scales to the SMBH and the other of rate rhc from direct infall of cold gas from galactic 
scales, TTiQs comes solely from direct infall of hot gas from galactic scales that remains in 
the hot component when reaching the SMBH. Suppose that the hot gas would have had 
an accretion of rate rhh feeding the SMBH, if it did not cool. But when the hot gas cools 
through its way to the SMBH scale from galactic scales, jfih will be split into rngs and rhhc- 
We are interested in computing rhqs here. 

First, a qualitative consideration. Assuming, conservatively, the no-outflow case where 
hot gas density scales with radius as p{r) oc r~^/^ and that the hot gas is at the virial 
temperature in the radial range r < r^, then, it can be shown that the ratio of gas cooling 
time to dynamical time is tcooi/tdyn T/A(T), where A(T) is the volumetric cooling 
function. Since A(T) oc T'^ with 7 < 1 (7 ~ 0.5 for T > lO^K with bremsstrahlung cooling) 
for the temperature range of interest here, tcooi/tdyn is an increasing function of decreasing 
radius from the SMBH in the radial range r < tb- On the other hand, in the radial range 
r > tb where rotation curve is nearly flat (Sofue & Rubin 2001; Gavazzi et al. 2007), 
assuming that infalling gas density goes as r~^, then, as long as 7 is greater than 1, the ratio 
tcooi/tdyn decreascs with decreasing radius, reaching a minimum at r^. As we have argued 
in Cen (2006), 7 is likely close to 2. The above two considerations indicate that tcooi/tdyn 
has a minimum at r^. Thus, one may estimate, probably to within a factor of a few, the 
amount of hot gas that will ultimately feed the SMBH, using physical arguments based on 
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the ratio tcooi/tdyn a-t tb- Guided by this, let us denote a critical hot QS-Bondi accretion 
rate as mg, at which tcooi/tdyn = 1 at r^. We use the cooling function from Sutherland & 
Dopita (1993), assuming solar metallicity for the hot gas. 

Let us now make a more quantitative calculation. There are two separate regimes in 
the hot accretion rate space. For m/j <^ mo, almost the entire hot gas accretion rate m/j 
may be expected to ultimately feed the SMBH in hot flow. In the other regime where rhh is 
not much smaller than mo, a progressively larger fraction of the hot gas will cool and "drop 
out" to fall onto the accretion disk to (potentially) feed the SMBH in cold flow or form 
stars. Under the assumptions (1) that the fraction of hot gas that "drops out" to become 
cold gas at each radius is directly proportional to the net cooling rate, (2) that already 
cooled gas does not get heated back up to return to the hot phase (i.e., assuming gas is 
thermally unstable, once being cooled) and (3) that the (remaining) hot gas is at the virial 
temperature at each radius, then, we may write the evolution equation for the hot gas: 

d{prr^Vr) 2fj,pf.r'^A(Tr) ApmprprCTrVr dlnTr pniprpvf dlnv^ 

dr ~ SkTrTUp SkTr dlnr S^T dlnr ' ^ ^ 

where r is the SMBH-centric radius, pr, Vr, (Tr and are the hot gas density, radial 
velocity, velocity dispersion and temperature at r, respectively. The first term on the 
right hand side of Equation (4) is the hot gas drop-out rate due to radiative cooling and 
the second term due to heating by the total gravitational energy released, with some of 
the released gravitational energy going to heat up the (remaining) hot gas to the virial 
temperature (the third term) and to accelerate the radial velocity (the fourth term). Let 
us define x = r/rs and y = {p/po){^'^/''^B)i'^r/vB), where po is the gas density at that 
corresponds to the critical accretion rate of mo, vb is the radial velocity at r^. Note that y 
is a "normalized" mass flux for the hot flow and would be constant in the absence of gas 
"drop-out". In Equation (4) we still need to specify Vr, which normally can be obtained by 
solving the radial velocity equation under the combined gravity and pressure force. But 
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here we will take a somewhat different approach by taking into account "non-spherical" 
components in our spherical Equation (4), to be consistent with our Assumption (3) above. 
That is, we require some viscous processes to turn radial kinetic energy into thermal 
energy, which, for example, may be accommodated by internal shocks in each radial shell. 
Hence, in conjunction with Assumption (3), we make one further simplifying Assumption 
(4) that Vr = (Jr, where Vr may be considered as an "assemble" radial velocity at r. In 
other words, some extra countering pressure force will be generated, if gas motion becomes 
supersonic. We note that the combination of Assumptions (3,4) precisely conserves energy 
in a Keplerian potential, which in our case corresponds to the radial range r < vb that 
primarily determines the final hot accretion rate (see below). However, it is noted that, if a 
different assumption, say, Vr — 0.5(7^, is made, the results only change ~ 20%. 

Thus, finally, we obtain a simplified hot gas evolution equation: 
dy ^ lAy^ HT,)TbvI 4y 4y dlnT, 

dx x2 K{Tb)T^vI ?,x Zx dhix ^ ' 

where [A(Ta;),A(TB)], (T^, Tb) and {vx, Vb) are the volumetric cooling functions, 
temperatures and radial velocities at (x, x — 1), respectively. In obtaining Equation (5) we 
have used the definition that tcooi/tdyn = 1 at x = 1 for the hot accretion rate of mo- The 
coefficient 1.4 refiects the difference between [p/nipY and ntUe that is used in association 
with A, where and ng are the total volume density of ions and electrons, respectively. 

Solving Equation (5) is a boundary condition problem. The outer boundary condition 
is provided by i/j, which corresponds to a hot gas density that would give a hot accretion 
rate of yjrho, if the hot gas is not allowed to cool. For a given value of we integrate 
Equation (5) inward from some large x to x — Q to obtain = y{x = 0). To facilitate a 
more tractable calculation, we assume, reasonably, that at r > the rotation curve is flat 
and at r < rs the rotation curve is Keplerian. Figure 1 shows y/ as a function of yi for five 
cases with Tb = (105•^ 10^ 10^ ^ 10^ lO^-^) K, corresponding roughly to M = 10^ - 10*^ Mq. 
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Fig. 1. — shows ?// as a function of initial i/i, computed based on Equation (5), for five cases 
withTs = (105•^ 10^ 10^•^ lO^ lO^-^), with (sohd, dotted, dashed, long-dashed, dotted-short- 
dashed) curves, respectively. 
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We find that, when > 1, = 0.6 ± 0.3, for Tg = lO^-^ - lO^-^K. The relatively weak 
dependence of rnqs^max on Tb (hence SMBH mass) reflects features in the cooling function 
A. At thh < O.lmo, practically all hot flow goes to the SMBH, as expected. 

It is prudent to note that our spherical treatment would break down at some small 
radius, because the angular momentum of the hot accreting gas is finite and the hot 
accretion may significantly deviate from sphericity at some small radii. It is thus instructive 
to check where the gas "drop-out" from the hot fiow practically stops. Figure 2 shows 
the evolution of Ux as a function of x with yi = 10^ (i.e., m/j = lO'^rho) for two cases with 
Tb — (10^-^, 10^-^) K. We see that most of the drop-out occurs at large radii x > 1. Then, 
most of the remaining drop-out within Bondi radius (i.e., x < 1) occurs at x > 0.1. This 
indicates that hot gas drop-out in the high i/i limit ceases at quite a large radius, close to 
the Bondi radius, where the assumption of a hot spherical flow may be relatively good. For 
lower values of i/i, this would occur at still larger radii. Therefore, our computed rate rriQs, 
while not exact, may be good to within a factor of order a few, which is adequate for our 
illustration. Figure 3 shows the computed maximum hot accretion rate onto the SMBH, 
rriQs^max, ^ a function of M, with the following simple fitting formula accurate to within 



shown as a dotted- long-dashed line in Figure 3. 

Our essential point is that, in general, the gas accretion onto SMBHs is composed 
of two concurrent accretion modes and in particular the hot accretion fiow of rate rriQs 
exists at all times. We note that the idea of hot accretion fiow around SMBHs is not 
new. Several well known classes of physical models have been proposed to understand 
the apparent radiative inefficiency of accretion onto some SMBHs at a low accretion rate: 
advection-dominated accretion fiow (ADAF; Ichimaru 1977; Narayan & Yi 1994; Kato, 



30% for the SMBH mass range M = 10^ - lO^-^ M©: 




(6) 
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Fig. 2. — shows function of initial x, for two cases with Tg = (lO^'^, 10'^, lO'''^) and 

both with an initial Ui = 10^ (i.e., rhh = lO^TfiQ^M)), with (solid, dashed) curves, respectively. 
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Fig. 3. — The solid and dashed curves show the maximum hot QS-Bondi accretion rate and 
mo, respectively, in Eddington units as a function of M. The dotted curve shows the Bondi 
accretion rate in a "dormant" state. Also shown as dotted-long-dashed curve is a simple 
fitting formula, Equation (6) 
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Fukue, & Mineshige 1998), convection-dominated accretion flow (CDAF; Quataert & 
Gruzinov 2000) and an inflow-outflow variant (Blandford & Begelman 1999), among others. 
A mere extension proposed here is that one should expect a co-existence of a hot accretion 
even in the regime where the total accretion rate is high (i.e., » rho). The concept of 
a possible co-existence of different accretion modes near the SMBH has been put forth 
by Chakrabarti & Titarchuk (1995), where the shock-heating-produced hot gas within 
ten Schwarzschild radii was proposed for the primary purpose of Comptonization of disk 
emitted photons. We self-consistently compute the accretion rate of the hot component in 
a way that takes into account the en-route processes. More importantly, here, we point 
out that the mix of cold and hot flow, in particular, the maximum hot flow rate rngs^max, 
depends strongly on M, which may have profound ramiflcations with regard to AGN jet 
production, among others. We stress that, while the numerical values obtained by solving 
the simplifled Equation (5) are certain to contain errors, the strong dependence of rnqs^max 
on M is primarily a result of the strong dependence of mo on M (dashed curve in Figure 3) 
that does not rely on all the assumptions employed to arrive at Equation (5). 



2.2. Hot Accretion Flow and AGN Jet Formation 

Inflow hot gas to SMBHs is likely signiflcantly pre-enriched with magnetic fluxes from 
earlier supernova remnants. Hence, realistic treatments of its dynamics would require 
explicit inclusion of magnetic fleld (e.g., Proga & Begelman 2003; Casse & Keppens 2004; 
Kato, Mineshige, & Shibata 2004; Proga 2005; De ViUiers et al. 2005; McKinney 2005; 
Krolik, Hawley, & Hirose 2005; Hawley & Krolik 2006). An emerging feature from these 
studies of inflow gas of various initial conflgurations with the inclusion of magnetic fleld is 
the strong possibility of the production of bipolar magnetic jets and funnel outflows. The 
former are akin to the magnetic tower (Lynden-Bell 1996,2003,2006). 
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One primary discovery that is borne out from these simulations is that, without 
significant external thermal pressure with effectively a large "scale-height", a magnetic 
tower fails to develop into a relativistic jet. The limited ability of magnetic self- 
confinement/collimation was previously known (e.g., Eichler 1993; Spruit, Foglizzo, & 
Stehle 1997; Okamoto 1999). We follow this indication and make the following critical 
ansatz: the QS-Bondi flow is responsible for the confinement and collimation of (twin) jets, 
by providing the requisite external pressure, and the strength (power) of the coUimated jets 
is directly commensurate with the external pressure hence the QS-Bondi accretion rate. 
This conjecture could find its precedent in the seminal work of Rees et al. (1982), where 
an ion-pressure supported thick disk was proposed to coUimate a pair of radio jets. We 
suggest here that the external pressure provided by the hot accretion can "naturally" act 
to coUimate and power the jet, as simulations have begun to indicate. Another important 
discovery from relativistic MHD jet simulations (e.g., McKinney 2005; Hawley & Krolik 
2006) is that the energy extraction efficiency, (3{a/M) (a is the spin of the SMBH), in the 
form of jet power in terms of accretion rate is a strong function of the spin of the central 
SMBH, in the vein of Blandford & Znajek (1977). Combining these ideas, we propose that 
the jet power goes as 

Pradio = (5{a/M)mQsMEddC^- (7) 

Recently, Allen et al. (2006) provide strong observational evidence that hot accretion 
may indeed be responsible for radio jet power in galaxies with moderate jet powers 
{Pjet ~ 10^^ — lO^erg/s). While our proposed picture was conceived before that paper 
was submitted, these new observations provide an impetus. Further, Wilson & Colbert 
(1995) ffist made the suggestion that SMBH spin may be responsible for the dichotomy 
between RQ and RL AGNs, although their overall scenario differs from ours (see below). 
McKinney (2005) provides a fitting formula p{a/M) = 0.07[nH/^Hia/M = where 
fin = etc/ {2MrH) is the rotation frequency of the BH and th = Tp(l + [1 — (a/M)^]^/^) is the 
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radius of the SMBH horizon for angular momentum J = aGM/c; quantitatively, P drops 
by a factor of 177 from a/M = 0.998 to a/M = 0.60 using this fitting formula. We assume 
a/M = 0.998 to be the maximal possible SMBH spin (Thorne 1974) and a/M = 0.60 
roughly corresponds to the SMBH spin without significant gas accretion (Volonteri et al. 
2005; V05 hereafter). Hawley & Krolik (2006) also provide a different fitting formula, 
j3{a/M) = 0.002(1 - a/M)'^. Note that the Bardeen-Petterson (1975) effect likely aligns 
the spin of the inner accretion gas with the spin of the SMBH. Thus, we consider the issue 
of inevitable mis-alignment between the large-scale gas angular momentum vector and the 
SMBH spin irrelevant to jet production efficiency at the benefit of ridding ourselves of one 
additional potential parameter, because the alignment time scale of the SMBH spin to large 
scale accretion disk is probably significantly shorter than the time scale for enough gas to 
be accreted via disk to reach a maximally spinning SMBH (e.g., Natarajan & Pringle 1998). 
A related consequence will be that the powerful radio jets are almost always perpendicular 
to the large-scale accretion disk or perhaps the galactic plane. 



2.3. Hot Accretion in "Dormant" SMBHs 

For the majority of time, a SMBH may be "dormant", i.e., inactive. Therefore, an 
estimate of the hot accretion in dormant SMBHs is useful. To do this, we need first to place 
the physical processes on SMBH scales in the context of galaxy formation, by beginning the 
discussion on galactic scales. In the standard cold dark matter model (Spergel et al. 2006) 
structures grow largely by mergers and accretion (e.g., Lacey & Cole 1993, LC hereafter). 
Keres et al. (2005), based on detailed analyses of smoothed particle hydrodynamic (SPH) 
simulations, have made an important discovery that there are two gas accretion modes 
on galactic scales, where massive galaxies in halos > 10^^ M© accrete predominantly hot 
diffuse gas, whereas lower mass galaxies accrete progressively more cold gas. Birnboim & 
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Dekel (2003; Dekel & Birnboim 2006) provided a plausible physical explanation for this 
phenomenon, proposing that in large halos the cooling time of gas in the post-shock region 
is longer than the gas compression time scale and hence shocked gas is unable to cool, 
resulting in hot gas. Given that, we would like to ask a somewhat different question: how 
does a massive halo maintain a hot quasi-spherical gaseous atmosphere? 

A gaseous atmosphere around a galaxy-size dark matter halo may be subject to thermal 
and Jeans instabilities. If we start with a just virialized halo with a hot atmosphere, the 
gas cloud is, by definition, (marginally) Jeans stable. Therefore, perturbations that are 
most susceptible to gravitational fragmentation are the largest possible ones of wavelength 
A = 2r, at radius r. We will assume that some heating processes effectively transmit 
energy to potentially counter cooling at a speed close to the sound speed, which will be the 
case if dissipation is through acoustic waves or low Mach number weak shocks. Then, the 
non- fragmentation condition is: tcooi/tsc > 1, where tcooi is the cooling time and tsc = 2r/cs 
the sound crossing time at r (c^ is the local sound speed). Figure 4 shows tcooi/tsc as a 
function of radius for halos of different masses. 

The existence of a critical halo mass. Ma ~ 10^^ M©, is evident from Figure 4. For 
halos more massive than Md, tcooi/tsc > 1 at all radii, indicating that gas at all radii is 
stable to fragmentation and the whole hot atmosphere remains quasi-spherical (at the 
virial temperature). As a result, subsequent inf ailing gas clouds may be shock heated and 
ram-pressure removed from their host dark matter halos, giving rise to "hot" accretion. 
But, even with the non-fragmentation condition met, some energy heating sources are still 
required in order to maintain the hot atmosphere in a quasi-steady state. For halos with 
M > Md, heating due to gravitational energy released by accreting satellite halos/galaxies 
may keep the gas roughly at the virial temperature, as shown by direct high resolution 
simulations (Naab et al. 2006). Additional heating, such as feedback from star and/or 
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Fig. 4. — shows the ratio of coohng time to sound crossing time as a function of radius for 
galaxies of different masses at redshift z = 6 (dashed curves) and z = (sohd curves); r„ 
is the virial radius. We use the NFW halo density profile (Navarro, Frenk, & White 1997) 
(top curves at small radii) with halo concentration parameter as a function of halo mass and 
redshift taken from Dolag et al. (2004) based on simulations, (using a different, Moore et al. 
(1998) density profile hardly reaches a different conclusion). We use the Komatsu and Seljak 
(2001) method to compute the gas density and temperature profiles, where the ratio of the 
gas density to total density at large radii is normalized to the global ratio. A metallicity of 
0.3 Zq is used for the halo gas and the cooling functions are taken from Sutherland & Dopita 
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SMBH formation, may also play an important role. In the absence of any heating source, 
the hot atmosphere will contract in the gravitational potential well. However, since the 
non- fragmentation condition is met, the contraction may be able to maintain the original 
density profile. 

For halos below there is a radial range where the gas cloud has cooling time that is 
sufficiently short to cause it to fragment and a quasi-steady coherent hot gas cloud at virial 
temperature at the fiducial gas density may not be maintained. Consequently, the amount 
of hot gas may be reduced until the remaining hot, more rarefied atmosphere satisfies the 
non-fragmentation condition. Gas accretion will thus progressively become more dominated 
by cold accretion in halos less massive than Ma due to a thinner atmosphere. Figure 5 
shows the relative reduction of hot gas density with respect to the fiducial case where the 
ratio of the gas density to total density at large radii is normalized to the global ratio; we 
have assumed for simplicity that the shape of the gas density profile remains unchanged. 
Taking the reduced hot gas density profile, as shown in Figure 5, we can estimate the 
QS-Bondi accretion rate onto the SMBH in dormant state, which is shown as the dotted 
curve in Figure 3. Note that this is only suggestive, since the overall atmosphere is subject 
to many other effects and likely dynamical. 

The results based on this very simple analysis are broadly consistent with simulations 
of Cox et al. (2004) and Keres et al. (2005), who find that gas infall is almost all cold at 
Mh < 10^^ M0, with the fraction of cold accretion dropping to ~ 50% at ~ 3 x 10" Mq, 
to nearly zero at Mh ^ 10^2 The assumed zero metallicity in Keres et al. (2005) may 
have underestimated Ma by a factor of about two (Cattaneo et al. 2006), which would 
render a somewhat still better agreement between simulations and our simple model. It 
should be stressed that, while our simple physical picture intends to capture the essence 
of the primary physics, it is clearly incapable of accounting for all complexities in real 
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Fig. 5. — shows the relative reduction of gas density with respect to the fiducial case as a 
function of halo mass, assuming that the shape of the gas density profile remains unchanged, 
for redshift z = (0, 2, 4, 6) with solid, dotted, dashed and long-dashed curves, respectively. 
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situations. For example, the merger/accretion processes are often anisotropic with a 
large fraction of activities (mass and energy flows) occurring along quasi-linear structures 
(filaments) that bridge dense lumps (Keres et al. 2005). But, for our purpose, it provides us 
with a simple quantitative way to represent the "outer boundary" and initial conditions for 
physical processes on smaller, SMBH scales, insofar as it fits direct simulations reasonably 
well. As a side note, it is interesting that remains little changed at z = 6 compared to 
at ^ = 0. Several competing factors contribute to this weak evolution while higher density 
at higher redshift gives shorter dynamical time [oc (1 + z)"^/'^] and still shorter coohng time 
[oc (1 + z)~^], lower concentrations of halos and a fixed halo mass corresponding to a higher 
virial temperature at higher redshift reduce the cooling time fortuitously by a compensatory 
factor. 

It is useful to clarify an apparent inconsistency here. A flat rotation curve in the radial 
range r > tb may be a very good approximation (Sofue & Rubin 2001; Gavazzi et al. 
2007). Thus, the gas density profile there may be relatively steep, possibly in the range r'^ 
to r~^, as the solution to Equation (5) in §2.1 seems to suggest. But here (Figure 4) we 
show that the same ratio increases with decreasing radius. The difference is that, in the 
"dormant" state, the gas density profile has a fiat profile at the center ("core"), while, in an 
"elevated" state, the central total gas density profile may be closer to an isothermal profile 
(e.g., Cen 2006), although the slope of the hot component becomes progressively shallower 
with increasing total accretion rate. 

2.4. What Triggers Powerful Radio Activities? 

We now consider how a dormant SMBH may be triggered to become active. We will 
advocate the mechanism proposed by Barnes & Hernquist (1991) & Mihos & Hernquist 
(1996) that significant galaxy mergers drive gas inward to fuel the central SMBHs, which in 
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some cases may trigger powerful jets, as we argue below. That powerful radio sources may 
be triggered by major mergers was pointed out long ago by Heckman et al. (1986), based 
on observational evidence. 

We will adopt some of the elements in V05 on how the evolution of SMBH spins is 
governed by mergers of SMBHs and gas accretion in the context of hierarchical structure 
formation. In the absence of gas accretion but in the presence of SMBH mergers, the spins 
of SMBHs roughly retain the value given at birth, which, as in V05, is assumed to be 
~ 0.60 and be the typical spin of each SMBH prior to a major merger. This means that 
elliptical galaxies, once formed, have just been eating up mostly smaller galaxies with their 
central SMBHs swallowing up smaller SMBHs randomly, which decreases a/M, if the initial 
spin is high. In other words, most elliptical galaxies get "reset" to a moderate spin during 
the long relatively "quiet" phases, if the initial spin is high, because they just enjoy a lot 
of minor mergers that spin it down but not significant gas accretion that would spin it 
up. As V05 have shown, major binary mergers of SMBHs with mass ratio q greater than 
0.125 tend to produce faster spinning SMBHs, while less major mergers typically spin down 
SMBHs. Gas accretion is very effective to spin up a SMBH: a maximally rotating state 
(a/M = 1) is reached when the final SMBH mass is equal to times the original SMBH 
mass (Bardeen 1970), i.e., roughly one e-folding, starting with a Schwarzschild SMBH. 
But, in the following, we intend to make some general statement without adhering to a 
specific model, such as in V05 where the amount of accreted gas for each major merger 
event is tightly pegged to the bulge rotation velocity (Volonteri, Haardt, & Madau 2003). 
In general, there may be four possible outcomes from significant mergers with respect to 
radio emission: 

(1) A relatively low-strength major merger does not produce maximally spinning 
SMBHs, even with both the binary SMBH coalescence and disk gas accretion, thus resulting 
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in a RQ AGN during its entire (limited) accretion history. While it is very uncertain, 
we tentatively denote mergers with q ~ 0.1 — 0.5 for this class and call them "significant 
mergers" . 

(2) The binary SMBHs coalesce relatively early in the major merger event (q > 0.5) 
and that significantly increases the spin of the remnant SMBH. Then, continuous gas 
accretion causes the SMBH to maximally spin at a certain time. Subsequently, the AGN 
becomes radio-loud. In this case, an initial RQ AGN phase may precede the RL phase. 

(3) A major merger {q > 0.5) drives disk gas accretion to gradually spin up the SMBH. 
Then, at a later time (when the disk gas accretion is still active) the binary SMBHs coalesce 
to cause the remnant SMBH to "switch" to the maximal spin. Subsequently, the AGN 
becomes radio-loud. In this case, an initial RQ AGN phase may also precede the RL phase. 

(4) A major merger drives gas accretion to spin up the SMBH until the maximum spin. 
But the binary SMBHs are unable to coalesce during the entire active phase. In this case, 
an initial RQ AGN phase may again precede the RL phase. The possibility of having three 
separate cases (2-4) is due to the still uncertain time it takes for two SMBHs to merge. But 
in all three cases, the triggers are truly "major mergers" . 

In our picture, RQ AGNs may be due to either a significant merger (Case 1) or the first 
phase of a major merger (Cases 2-4). But, since the frequency of merger events decreases 
sharply with increasing merger strength (i.e., mass ratio) (e.g., LC), we expect that RQ 
AGNs are likely to be predominantly contributed by Case (1). On the other hand, RL 
AGNs may be predominantly in the latter stage of the AGN event when the SMBH mass 
has roughly increased by one e-folding (Cases 2-4). 
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3. Some Observational Tests and Predictions 

3.1. Strong Dependence of McLximum Jet Power on Black Hole Mass 

Equation (7) and Figure 3 combine to yield the maximum strength of radio jets that 
is a strong function of M, shown in Figure 6 for a Kerr BH with a/M = 0.998 (red 
curves) and a/M = 0.60 (blue curves), respectively, using the fitting formula for (3 {a/M) 
from McKinney (2005). Because ttiqs is expected to be non-zero for any SMBH, even in 
"dormant" state (the dotted and long-dashed curves in Figure 5), radio jets should emerge 
near SMBHs of all galaxies, as supported by recent observations (Ho 1999; Ulvcstad et al. 
1999; Kukula et al. 1999; Ulvestad & Ho 2001; Nagar, Wilson, & Falcke 2001; Giovannini 
et al. 2001; Bicknell 2002; Mundell et al. 2003; Falcke, Kording, & Nagar 2004; Anderson, 
Ulvestad, & Ho 2004; GaUimore & Beswick 2004; Nagar, Falcke, & Wilson 2005; Filho, 
Barthel, & Ho 2006; GaUimore et al. 2006). 

However, the strength (and probably scale) of radio jets varies vastly, seen in Figure 
5. For example, for low mass SMBHs (< 10^ M©), the maximum radio power is predicted 
to be low, in the range of ~ 10^^ — lO^^erg/s, in agreement with observations (e.g., Greene, 
Ho & Ulvestad 2006). Seyfert galaxies, having low M SMBHs, should not display powerful 
jets, as known (e.g., Lai, Shastri, & Gabuzda 2004). At the high end M > 10®-^ - 10^° Mq, 
radio jets could reach a maximum power in the range of 10^^'^ — 10^^'^erg/s for a nearly 
maximally rotating Kerr BH, fully consistent with observed radio jets (Fanaroff & Riley 
1974; Laing, Riley, & Longair 1983; Wall & Peacock 1985; Spinrad et al. 1985; McCarthy 
1995). The maximum radio power roughly scales as (indicated by the dot- long-dashed 
line in Figure 6), a very steep dependence on M that is in fact in good agreement with 
observations (Ledlow & Owen 1996; Laor 2000; Lacy et al. 2001; Boroson 2002; Jarvis 
& McLure 2002, McLure & Jarvis 2004; Metcalf & Magliochetti 2005; Best et al. 2005). 
We find that, with the adopted /9(0.998), the following formula provides a good fit for the 
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Fig. 6. — shows the maximum radio jet power as a function of M for SMBHs for two cases 
with a/M = 0.998 (sohd curve) and a/M = 0.50 (dashed curve), respectively. Also shown as 
dotted and long dashed curves are the counterparts with a fiducial gas density profile, i.e., in 
a "dormant" state. The dot-long-dashed curve indicates Pradio = 10^'^'^(M/10^ M0)^erg/s. 
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maximum jet power as a function of M: 

shown as dot-long-dashed hne in Figure 6. However, given the uncertainties with respect to 
both (3 and rnqs^maxi Equation (8) might be accurate only to an order of magnitude. 

Finally, we show the corresponding radio power from SMBHs in dormant state in Figure 
6 as the dotted and long-dashed curves, respectively, for a/M — 0.998 and a/M — 0.60. We 
see that, even in "dormant" state, a substantial radio power (> lO^^erg/s) from the nucleus 
should be expected for SMBHs with M > IO^Mq, while for low M < IO^Mq SMBHs the 
radio power is expected to be smaller than ~ lO^^erg/s. 



3.2. Evolution of Radio-Loud AGN Fraction with Redshift 

Our scenario suggests that the RQ AGNs may be predominantly produced by 
significant mergers that have lower strengths than major mergers that result in RL AGNs. 
This property may be manifested in the evolution of the relative fractions of the RL and RQ 
AGNs, because the distribution of merger strength with redshift may be precisely predicted 
in the CDM model and is expected to evolve. We estimate this using the Extended 
Press-Schechter (EPS) theory (Bond et al. 1991; Bower 1991; LC). We use the instantaneous 
merger rate for a halo of mass Mp to merger with a halo of mass AM = Mf — Mp to form a 
halo of mass Mf (LC): 



R.iMp ^ Mf)dMf ^ ^^(TT^iy^exp 



dMf, (9) 



where u/ and are the density variances of a sphere in the linear density field extrapolated 
to 2; = that contains mass Mf and Mp, respectively; 6c{t) is the critical density for collapse 
(equation 2.1 of LC); t is age of the universe at the epoch examined. 
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Fig. 7. — shows the ratio of merger rate with mass ratio g > 0.5 to those with g > 0.1 
for a halo mass of 5 x 10^^ Mq. The standard WMAP3-normahzed cold dark matter model 
(Spergel et al. 2006) is used: Qm = 0.26, A = 0.74, Hq = 72km/s/Mpc, n, = 0.95 and 
as = 0.77. 
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In the absence of more quantitative understanding of the physics of gas accretion 
and radio emission and their dependences on merger strength and SMBH mass, among 
others, it is uncertain to compute the dependence of RL fraction (RLF; Jiang et al. 2006) 
of quasars on the mass/luminosity of galaxies. For this reason, we will only compute the 
redshift evolution of the relative rates of mergers of different strengths at a fixed halo 
mass, hoping that there is little evolution in accretion processes at a fixed SMBH mass and 
merger strength. Figure 7 shows the ratio of merger rate with mass ratio q > 0.5 to those 
with g > 0.1 for a halo mass of 5 x 10^^ Mq. If we use this ratio as a proxy for the RLF, we 
find that RLF of quasars is a decreasing function of increasing redshift, a trend that is in 
agreement with observations (Schmidt et al. 1995; Jiang et al. 2006). A good quantitative 
agreement between our model and observations may be obtainable, if one fine-tunes the 
mass ratio cutoffs for significant and major mergers. 

3.3. Is There A Dichotomy in AGN Radio Loudness Distribution? 

In our proposed picture there are two separate populations of AGNs. The RQ 
population that is primarily created by significant mergers likely substantially outnumbers 
the RL population that is triggered by more rare, major mergers. The peaks of the two 
populations may be well separated in radio loudness by, to the zero-th order, the ratio of 
(3{loud)/ (3{quiet), which may amount to 100 — 1000, if we identify Piloud) with the case 
with maximal rotation and (5{quiet) with the case with a/M ~ 0.5 — 0.6. There are, of 
course, a likely spread in a/M distributions. This is quite uncertain presently. 

Thus, both scatters in optical luminosity at a fixed radio luminosity and in radio 
luminosities at a fixed optical luminosity would likely help smear out both peaks to some 
extent in the radio-loudness axis. More importantly, the sum of two decreasing monotonic 
functions of radio-loudness each spanning a range that is larger than the displacement of the 
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two functions in radio-loudness will likely largely produce a smooth overall distribution that 
suppresses or completely removes the secondary peak at the high radio-loudness. However, 
if one restricted oneself to an AGN sample where the range of SMBH mass is comparable to 
or narrower than the displacement between the RL and RQ population, then the secondary, 
smaller peak at the high radio-loudness may become visible. 

We will use a simple toy model to illustrate this. Since the maximum radio jet power 
roughly scales as M^, then, if we assume that the optical luminosity scales as M (i.e., 
at a fixed disk accretion rate in Eddington units), the radio-loudness parameter R may 
be linearly proportional to M, albeit with a significant dispersion. Thus, we will adopt 
a standard Schechter function for the underlying radio-loudness distribution of the RQ 
population: nRQ{R)dR = A{R/ R^)^^-"^ eicg{—R/ R^)dR. The RL population is then assumed 
to bear the form: url^R) = Q.lnRQ^R/rj), where rj is taken to be 200 /3(0.998)//3(0.6)]. 
Observatioally, a cutoff in some apparent optical magnitude applied to a flux-limited 
sample essentially translates to a lower cutoff in M to both the RQ and RL populations. 
The corresponding lower radio-loudness cutoffs to the RQ and RL populations will be Rcut 
and rjRcut, respectively. But the exact form at lower cutoff in R is uncertain and we will 
use simple functional forms exp(— (i?cut/-R)''^^^) and exp(— (jyi^cutZ-R)^'^^) for the RQ and 
RL populations, respectively, which multiply their respective n{R) distributions. Figure 8 
shows four cases with different Rcut- While we cannot ascertain the results precisely, the 
trend appears to be consistent with the observed trend found in White et al. (2007) in that 
a small secondary peak is visible in their full sample, which is then signiflcantly enhanced 
when a high-redshift sub-sample is used. The reason is that taking only the high-redshift 
sub-sample may be effectively imposing a higher lower-limit on the SMBH mass. Similarly, 
the more prominent secondary peak found by Ivezic et al. (2002) may be a result of applying 
an optical magnitude cut {i* < 18.5) to their SDSS sample of AGNs. 




Fig. 8. — shows the radio loudness distributions for four cases with low cutoff equal to (1) 
no cutoff (sohd), (2) O.li?* (dotted), (3) i?* (dashed) and (4) lOi?* (long-dashed). 
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Therefore, our conclusion is that, while there are truly two separate populations of 
AGNs, one RQ and the other RL with the overall shift in the radio-loudness possibly 
amounting to a factor of 100 — 1000, a perfect AGN sample with an infinitely high 
observational sensitivity would in fact result in a radio-loudness distribution that is not 
bimodal. However, in an observational sample of AGNs with significant observational 
selection effects, primarily those which remove fainter optical objects, one may be able to 
unravel, incidentally, the two underlying populations in the form of a secondary peak at high 
radio-loudness. A definitive test of this picture may be obtained, for example, by applying 
the analysis of Ivezic et al. (2002) to a fainter optical cutoff, where the secondary peak 
should progressively shift to lower radio-loudness with a diminishing relative amplitude. 



3.4. Some Additional Properties and Tests 

Numerous other consequences and predictions may be made based on this general 
model. Here, we highlight several of them. 

First, both RQ and RL quasars are expected to be drawn from the same underlying 
elliptical galaxy population, in agreement with observations (e.g., Bettoni et al. 2001; 
Scarpa & Urry 2001; Kauffmann et al. 2003; Urry 2004). However, the average SMBH mass 
of RL quasars may be higher than that of RQ quasars by roughly one e-folding, consistent 
with observations (e.g., Corbin 1997; Laor 2000; Lacy et al. 2001; Boroson 2002; Jarvis & 
McLure 2002; Metcalf & Maghochetti 2005). 

Second, since the hot accretion and cold disk accretion are practically "decoupled" 
in the sense that varying the hot accretion may not significantly alter the disk accretion 
physics, RQ and RL quasars should have nearly identical average IR-optical-UV properties 
which are assumed to be primarily due to the disk accretion, in agreement with observations 
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(e.g., Francis et al. 1993; Zheng et al. 1997). 

Third, if each merger event that triggers an AGN phase is accompanied by a (central) 
starburst (e.g., , Barnes & Hernquist 1991), the SN la rate may be higher in RL AGNs 
than in RQ AGNs. This is simply because RQ AGNs are formed by less major mergers 
(Case 1 in §2.4) and therefore they will experience less star formation hence fewer SN la's. 
For Cases 2-4 the RL AGN phase may be delayed by roughly an e-folding or Salpeter time 
compared to the earlier RQ AGN phase. This delay time is shorter than the typical lifetime 
of progenitors of SN la supernovae of ~ IGyr and hence the SN la rate is rising during the 
period, which would also result in RL AGN having higher SN la rate. Both scenarios thus 
yield a trend that appears to be supported by observations (e.g, Delia Valle et al. 2005). 

Fourth, under the assumption that more major mergers produce more pronounced 
central stellar "cores" in the final eUiptical galaxy, carved out by the coalescence of the two 
central SMBHs (Faber et al. 1997; Milosavljevic et al. 2002), RL AGNs may preferentially 
be found in core elliptical galaxies. There is some observational indication that this may 
indeed be the case (de Ruiter et al. 2005; Capetti & Balmaverde 2006). If the other finding 
of Capetti & Balmaverde (2006), that RQ AGNs are found only in power-law galaxies, holds 
up, this suggests that RQ AGNs arc primarily produced by significant mergers (Case 1 in 
§2.4), as expected in the standard cold dark matter model because mergers rate generally 
decreases rapidly with merger strength (mass ratio). The intriguing link between RL 
quasars and core elliptical galaxies may be indicative that massive binary SMBHs merge on 
a relatively short time scale, perhaps shorter than lO^yrs (e.g., Zier and Biermann 2001; Yu 
2003; Escala et al. 2004). 

Fifth, since the radio power is a strong function of M {Pjet oc M^), among radio- 
powerful AGNs there should be correlation between clustering strength and radio power for 
radio-powerful AGNs, in agreement with observations (e.g., Overzier et al. 2003), under the 
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assumption that more massive galaxies are more strongly clustered. However, some AGNs 
with very massive SMBHs may be RQ. As a result, for powerful AGNs (i.e., mostly above 
some high M threshold) radio-loud and RQ samples may exhibit comparable clustering 
strengths and inhabit similar environments, as indicated by observations (McLure & Dunlop 

2001) . On the other hand, spatial clustering of a random set of RQ AGNs would be an 
"average" over all SMBH masses. As such, it is "diluted" and should in general be less 
strongly clustered than that of powerful RL AGNs, in agreement with observations (e.g.. 
Cress et al. 1996; Hall et al. 2001; Maghocchetti et al. 2004; Brand et al. 2005; Zheng et al. 
2006; Bornancini et al. 2006). 

Sixth, for Cases (2-4) (§2.4) it is possible sometimes that the cold disk accretion 
becomes diminished at some later stage, when the spin of the SMBH has become maximal 
and the hot accretion still remains vigorous. This may result in radio galaxies that do not 
have a significant amount of disk material, presumably including the purported obscuring 
torus. There are some observational indications for the existence of this type of objects 
(Whysong & Antonucci 2004; MuUer et al. 2004; Kharb & Shastri 2004). 

Seventh, in all Cases (1-4), at some later time the disk accretion may finally taper 
off. However, the hot accretion would gain in relative strength when total accretion 
or disk accretion decreases. By this time, on the order of a Gyr has elapsed since the 
onset of the merger event. We propose that this "traihng" AGN phase may be identified 
with Low-luminosity AGN (LLAGN) and LINERs (Low Ionization Narrow Emission-line 
Regions; Heckman 1980; Ho, Filippenko, & Sargent 2003). One specific prediction is that 
these objects are hkely to be radio- loud, supported by observations (e.g.. Ho 1999). More 
generally, there should be an anti-correlation between radio-loudness and disk accretion 
rate, in good agreement with observations (e.g., Oshlack, Webster, & Whiting 2001; Ho 

2002) . In addition, there may be a correlation between star formation activity and strength 
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of the AGN activities in these objects, since the later is on a dechning slope because of 
declining accretion rate and the former declines with time because of passive aging. 

Eighth, since major mergers, while rare, could occur multiple times in a Hubble time 
for some galaxies, restarting of radio jets may thus be expected, and observations support 
this expectation (e.g., Lara et al. 1999; Marecki et al. 2006). Because accretion disks from 
multiple mergers are generally not expected to be aligned, complex radio galaxy structures 
may be produced (e.g., Hogbom & Carlson 1974; Ekers et al. 1978; Leahy & Parma 1992; 
Merritt & Ekers 2002; Dennett-Thorpe et al. 2002; Gopal-Krishna, Biermann, & Wiita 
2003). 

3.5. A Thought on Radio Jet Feedback Effect 

The problem of the jet feedback is unsolved. However, if we take the cue from the sizes 
of powerful radio jets (PR Is and PR lis) (e.g., Urry & Padovani 1995), which often reach 
lOOkpc-lMpc, we may find some hint. At lOOkpc for a medium of sound speed of 300km/s, 
the sound crossing time is ~ 3 x lO^yrs. If most of the jet energy channels to and dissipates 
on these large-scales features, it means that the feedback of the jets on the surrounding gas 
will only be felt after this time scale. While this argument may be quite reasonable for 
powerful radio jets, one needs to be careful that the exact radial energy dissipation profile 
of jets on the surrounding medium is unknown. Purthermore, one might reason that less 
powerful jets in smaller SMBHs may start dissipating at smaller radii. But the countering 
effect that the surrounding gas medium will be proportionally less dense in less powerful 
jets, recalling that Pradio oc tHqs oc Phot(^)- As a result, it is unclear that jet feedback effect 
on less powerful jets may be more profound. As perhaps supportive evidence that the jets 
will be "slowed" down and dissipate by a denser medium at a more rapid rate is that the 
observed relic radio cavities in the cores of clusters of galaxies, where gas density is higher 
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than for example in regions surrounding isolated galaxies, are observed on a smaller scale 
of tens of kpc instead (e.g., McNamara et al. 2000,2001; Fabian et al. 2000; Blanton et al. 
2001). We will come back to examine this issue in much greater detail subsequently. 

4. Conclusions 

We suggest a co-existence of hot and cold accretion flows for any SMBH and propose 
that the hot accretion is a primary determinant of the radio emission from AGNs. Physically 
this is based on the notion that the external pressure provided by the hot accretion is 
required for coUimating the radio jets. The spin of the SMBH also plays a major role in 
terms of the efficiency by which energy may be extracted in the form of radio jets from hot 
accretion. We propose that the jet power is directly proportional to the hot accretion rate 
and spin-dependent energy extraction efficiency in the form of jets (Equation 7). 

A variety of predictions may be made and all seem to be in agreement with extant 
observations. Specifically, 

(1) Radio jets should emerge in all SMBH to varying strengths. 

(2) The maximum radio jet power increases with increasing SMBH mass approximately 
as M^. 

(3) While there two two separate underlying populations of RQ and RL AGNs, the 
apparently observed two peaks in the AGN radio loudness distribution are due to selection 
effect, caused by an imposed optical magnitude limit. 

(4) The RL fraction of quasars is expected to decrease with redshift in the cold dark 
matter model. 

(5) The host elliptical galaxies of RQ and RL quasars are drawn from the same 
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underlying elliptical galaxy population, although RL quasars are hkely to possess SMBH 
that are somewhat more massive than their RQ counterparts and may reside predominantly 
in core elliptical galaxies. 

(6) Low-luminosity AGN and LINERs may be the long, declining "trailing" phase 
following the initial, more luminous AGN phase. They should be RL and conform to the 
general, broad anti-correlation between radio-loudness and disk accretion rate. A correlation 
between star formation activity and strength of the AGN activities in these objects may be 
expected. 

(7) RL galaxies may be expected to be more abundant in Type la supernovae. 

(8) Among RL AGNs a correlation between radio power and clustering strength is 
predicted. 

(9) Radio-quiet and RL AGNs should have nearly identical average IR-optical-UV 
properties. 
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